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Abstract 

We report results from Monte Carlo simulations of systems of magnetic dipoles that relax through quantum 
tunneling, much as Feg crystals at very low temperature. For short times, a hole develops in suitably defined 
magnetic field pseudo-distributions, which matches the shape of the tunnel window (TW). Much later, ordinary 
field distributions P{h) develop similar holes if thermal energies are not much larger than the TW's energy. Still 
later, below the long-range ordering (LRO) temperature, P{h) exhibits the signature of LRO. 
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Magnetic relaxation in crystals of single- 
molecule magnets, such as Feg, has recently be- 
come a subject of great interest. Experimentally 
observed relaxation at temperatures T that are 
below some 10~^ of single-ion crystal anisotropy 
barrier energies U is temperature independent, 
and is duly attributted to magnetic quantum 
tunncUng (MQT) [1]. As explained in Ref. [2], 
hyperfine interactions with nuclear spins open up 
tunneling windows (TW's) of energy e^u that are 
comparable to magnetic dipolar energies. Then, 




(1) 



is the tunneling rate for spins at very low tem- 
perature, where T is some rate, r](x) ~ 1 for | 
X 1^ 1; ~ for a; > 1, and 2eh is the Zee- 
man energy change upon tunneling. From here on, 
times are given in terms of F"-'^. The results re- 
ported here in Figs. 1, 2(a), and 3 are for ri{x) = 1 
for I a; |< 1 and r]{x) — for | x |> 1; in Fig. 
2(b) we show results for ri{x) = exp(— | a: |). We 
shall refer to the former TW as flat and to the 
latter one as exponential. An unusual magnetiza- 
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Fig. 1. f /h versus h for the shown times after quenching 
and applying a field H = 2. All data points are for averages 
over 4 X 10^ independent runs for = 0.1. In the inset, • 
and o are for P_i_ and P_. 



tion M relaxation of fully polarized Feg like sys- 
tems, Mt=o — Mt (X -\/t, was predicted in Ref. [2] 
and later observed [3]. In an important piece of 
work, Wernsdorfer et al. [4] found that Mt oc ^/t 
if a small magnetic field H is applied at i = 0, af- 
ter quenching from a higher temperature. This be- 
havior has been explained recently [5] . In addition, 
since dMt{H)/dt = -2T J dhf{h,t)r]{h + H) at all 
times, where f{h,t) = [P_|_(/i,t) — P^{h,t)], and 
P+{h) [P-{h)] is the density distribution for fields 
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h acting on up [down] spins, /(/i, t) can be deter- 
mined by a suitable application of magnetic fields. 
[/ itself is not normalized, and is therefore not a 
distribution. In fact, it turns out that f{h) is pro- 
portional to the system's energy [5]]. A "hole" of 
"intrinsic width" , was thus observed to develop in 
time, and conjectured to correspond to the width 
in r]{h), i.e., to the width of the TW [4]. 

We report Monte Carlo simulations of Ising sys- 
tems of interacting dipoles on SC lattices. Initially, 
spins are allowed to flip readily at temperature 
Ta > C//10 > To (To ~ 2.5 is the LRO tempera- 
ture) until the energy per spin reaches some value 
Sa- A field H is applied then and spins only flip 
at rate F' [see Eq. (1)] thereafter, thus simulating 
temperatures below some lO^^U. We have previ- 
ously found [5] that Mt oc\ Sa \ t. Here, we re- 
port results for the evolution of the field pseudo- 
distribution function f{h) and ioi p{H) = + 
P-{h). Fields and energies are given in terms of 
nearest neighbor dipolar values. 

How P+{h) and P-{h) split, leading to a non- 
vanishing f{h), after thermalization at high T is il- 
lustrated in the inset of Fig. 1. The hole that devel- 
ops in f{h) after quenching and applying field H is 
shown in Fig. 1 . The shape of the hole is shown in 
Figs. 2(a) and (b), for flat and exponential TW's, 
respectively, for H = 2. Clearly, not only is the 
hole's width the same as the width of rj, as conjec- 
tured in Ref. [4], but f{h)/h matches ri{h) for the 
three shortest times shown in Figs. 2(a) and (b). 
When t >1, the hole widens, as shown in Fig. 2(a) 
for t = 4. 

On the time scale on which the hole in ///i de- 
velops, P{H) remains approximately independent 
of time. On the other hand, a hole which resembles 
the TW grows at much later times in P{H) if T < 
£w [6]. This is exhibited in Fig. 3. Still later, below 
the long-range ordering (LRO) temperature, P{h) 
exhibits the signature of LRO [7] whether T < 
or not. If, on the other hand, the energy of the sys- 
tem is kept constant, no LRO ensues, as expected. 
This is illustrated in Fig. 3. 
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and < are for t = 0, 0.0625, 0.25, 0.5 and 4, respectively. 
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Fig. 3. The field distribution P versus h for the shown 
times, and T = 0.1 or constant energy. In all cases 
Sa = -0.58, = 0.2 and H = 1. 



2 



